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ABSTRACT - The effect of heat source/sink on steady two-dimensional laminar MHD mixed convection 
boundary layer flow past a vertical permeable surface embedded in a porous medium saturated by a nanofluid is 
performed in this study. Numerical solutions of the similarity equations are obtained using the shooting method. 
Three types of metallic or nonmetallic nanoparticles, namely Copper (Cu), Alumina (Al 2 Os) and Titania (TiOf) 
are considered by using a water-based fluid to investigate the effect of the solid volume fraction or nanoparticle 
volume fraction parameter U of the nanofluid. The numerical results of the skin friction coefficient and the 
velocity profdes are presented and discussed. It is found that the imposition of suction is to increase the velocity 
profiles and to delay the separation of boundary layer, while the injection parameter decreases the velocity 
profiles. On the other hand, the range of solutions for the injection case is largest for Al 2 0 3 nanoparticles and 
smallest for Cu nanoparticles. 

Keywords - Boundary Layer, Nanofluid, Mixed Convection, Permeable Surface, Porous Medium, Magnetic 
field 



I. INTRODUCTION 

The research of magnetohydrodynamic (MHD) incompressible viscous flow has many important 
engineering applications in devices such as power generator, the cooling of reactors, the design of heat 
exchangers and MHD accelerators. Yih (1997) considered the effect of uniform blowing/suction on forced 
convective heat transfer of magnetohydrodynamic Hiemenz flow through porous media invoking the model of 
the porous medium proposed by Raptis and Takhar (1987). 

The problem of convection in a porous medium provides one of the basic scenarios for heat transfer 
theory and thus is of considerable theoretical and practical interest and has been extensively studied. Excellent 
reviews of the topic can be found in the books by Nield and Bejan (2006); Pop and Ingham (2001); Ingham and 
Pop (2005) and Vadasz (2008). The most basic problem for natural or free convection in a porous medium past 
a vertical flat plat was first studied by Cheng and Minkowycz (1997). There are several numerical studies on the 
mixed convection in a porous media and we mention here those by Harris et al. (2009); Rosali et al. (2011); 
Imran et al. (2012) and Mukhopadhyay (2012). On the other hand, nanofluids are engineered by suspending 
nanoparticles with average size below 100 nm in traditional heat transfer fluids such as water, oil and ethylene 
glycol. Fluids such as water, oil and ethylene glycol are poor heat transfer fluids, since the thermal conductivity 
of these fluids play important role on the heat transfer coefficient between the heat transfer medium and the heat 
transfer surface. Choi and Guarino (1995) showed that the addition of small amount (less than 1% by volume) of 
nanoparticles to conventional heat transfer liquids increased the thermal conductivity of the fluids up to 
approximately two times. Therefore, the effective thermal conductivity of nanofluids is expected to enhance 
heat transfer compared to the conventional heat transfer liquids. Some numerical and experimental studies on 
the forced and natural convection using nanofluids related with differentially heated enclosures have been 
considered by Jou and Tzeng (2006); Tiwari and Das (2007); Abu-Nada (2008); Oztop and Abu-Nada (2008); 
The temperature-dependent heat generation/absorption in energy equation is considered by Vajravelu and 
Rollins (1992). They adopted the model of Foraboschi and Federico (1964) to study the heat transfer in an 
electrically conducting fluid over a stretching surface. 

In many problems, there may be an appreciable temperature difference between the surface and the 
ambient fluid. This necessitates the consideration of temperature-dependent heat sources or sinks which may 
exert strong influence on the heat transfer characteristics. The study of heat generation or absorption in moving 
fluids is important in view of several physical problems such as fluids undergoing exothermic or endothermic 
chemical reactions. Exact modelling of the internal heat generation or absorption is quite difficult, some simple 
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mathematical models can express its average behaviour for most physical situations. Heat generation or 
absorption has been assumed to be constant, space-dependent or temperature dependent. Sparrow and Cess 
(1961) considered temperature dependent heat absorption in their work on steady stagnation point flow and heat 
transfer. Watanabe(1991) investigated the forced and free mixed convection boundary layer flow with uniform 
suction or injection on a vertical flat plate. Vajravelu and Nayfeh(1998) reported on the hydromagnetic 
convection at a cone and a wedge in the presence of temperature-dependent heat generation or absorption 
effects. Chamkha(1992) reported the non-darcy fully developed mixed convection in a porous medium channel 
with heat generation/absorption. Chamkha et al. (2003) discussed effects of thermal radiation on MHD forced 
convection flow adjacent to a non-isothermal wedge in the presence of a heat source or sink. Azim et al. (2010) 
reported viscous Joule heating MHD-conjugate heat transfer for a vertical flat plate in the presence of heat 
generation. Alam et al. (2007) investigated the similarity Solutions for hydromagnetic free convective heat and 
mass transfer flow along a semi-infinite permeable inclined flat plate with heat generation and thermophoresis. 
Chamkha et al. (2004) studied the mixed convection flow over a vertical plate with localized heating and 
cooling along with magnetic effects and suction/injection effects. Kandasamy et al. (2010) discussed the effect 
of temperature-dependent fluid viscosity with thermophoresis and chemical reaction on MHD free convective 
heat and mass transfer over a porous stretching surface in the presence of heat source/sink. Recently, Sharma et 
al. (2010) investigated the combined effect of magnetic field and heat absorption on unsteady free convection 
and heat transfer flow in a micropolar fluid past a semi -infinite moving plate with viscous dissipation using 
element free Galerkin method. 

In the present paper, we have considered the heat source/sink effect on MHD mixed convection 
boundary layer flow saturated by electrically conducting nanofluids over a vertical permeable flat plate with 
linear wall temperature proportional to the distance in a porous medium. 



In this paper the effect of heat source/sink on the steady two-dimensional laminar MHD mixed 
convection boundary layer flow past a vertical semi -infinite plate embedded in a porous medium filled with a 
nanofluid is studied. It is assumed that the free stream velocity and the ambient temperature (far flow from the 

plate) are U x andr x , respectively. It is also assumed that the temperature of the plate isr , where 
T > T r corresponds to a heated plate (assisting flow) and T < T r corresponds to a cooled plate (opposing 

flow). It is also assumed that the convecting fluid and the porous medium are in local thermodynamic 
equilibrium, the viscous dissipation is neglected, the physical properties of the fluid except the density are 
constant and that the Boussinesq approximation holds. Following the nanofluid equation model proposed by 
Tiwari and Das (2007) along with the Boussinesq and boundary layer approximations, it is easy to show that the 
steady boundary layer equations of the present problem are (Ahmad and Pop, 2010): 



II. 



PROBLEM FORMULATION 
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dx dy 
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(2) 



(3) 



Subject to the boundary conditions: 



a=0, p = tfc ,T = T w 0t) at y = 0 



J-p 7„ as y-f en 



(4) 
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Integrating Equation 2 with the boundary conditions Equation 4, it becomes: 

*t H *T ' P (5) 

Here, x and y are the Cartesian coordinates measured along the plate and normal to it, respectively, u 
and v are the velocity components along x and y axes, respectively, T is the temperature of the nanofluids, g is 

the acceleration due to gravity, v ( x) is the mass transfer velocity with v ( x ) < 0 for suction and 
v w ( x) > 0 for injection, <j> is the nanoparticle volume fraction, p . is the dynamic viscosity of the base 
fluid, p f and J3 are the coefficients of thermal expansion of the fluid and of the solid, respectively, p and 
p are the densities of the fluid and of the solid fractions, respectively, p . is the viscosity of the nanofluid 
and a , is the thermal diffusivity of the nanofluid; k is the permeability of the porous medium; a is the 

electrical conductivity; B„is the externally imposed magnetic field in the y-direction. The induced magnetic 

field effect is neglected for the small magnetic Reynolds number flow, by Shercliff (1965). It is also assumed 
that the external electric field is zero, the electric field owing to polarization of charges and the Hall effect are 

neglected; Q is the volumetric rate of heat generation/absorption; C is the heat capacity at constant pressure, 

and the Thermophysical properties of fluid and nanoparticles are given by Oztop and Abu-Nada (2008): 

- - K *l 

^ ~ a _ ^ ■ a v--( P c v ) n} 

(6) 

% _ (K s +2K f )-2tp(K f -K 3 } 
K f ~ (K s + 2K / ) + <p(K f -K s ) 

Where, q ) is the heat capacity of the nanofluid as expressed by Khanafer et al. (2003) and Abu-Nada 
(2008). The viscosity of the nanofluid p. can be approximated as the viscosity of a base fluid p f containing 

dilute suspension of fine spherical particles. The effective thermal conductivity of the nanofluid k. is 

approximated by the Maxwell-Garnett's model, which is found to be appropriate for studying the heat transfer 
enhancement using nanofluids (Abu-Nada, 2008). Similarity solutions of Equation 3 and 5 subject to the 
boundary conditions (4) of the following form: 



1/2 

(7) 



Where, Pe x = U x x I a f is the local Peclet number for the porous medium and y is the stream function, 
which is defined in the usual way as u = dy/ / 8y and v = 8 y/ fdx . Thus, we have Equation (8): 



(8) 
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Where, primes denote differentiation with respect to 77 . In order that Equation 1 to 3 subject to the boundary 
conditions (4) admits a similarity solution, we have to consider that v (x)has the following expression 
Equation (9): 



la f (2Pe,) 1/z 



(9) 



where, / is the constant mass transfer parameter with / > 0 for suction and f 0 <0 for injection. 

Substituting Equation 6 and 7 into Equation 3 and 5, we obtain the following system of ordinary differential 
equations: 



r = 1 + 



p/Ap> 



W — M z {f' — 1) 



(10) 



Along with the boundary conditions: 

f(0) = f 0 1 0(0) = 1, /'(«)-» 1 , 0(«O -» 0 

Where, A is the constant mixed convection parameter which is defined as Equation (12): 

- HfL. 

A = — - 



(11) 



(12) 



With Ra t = p f gk p f (T w ~ T r )x j 'a f ju f being the local Rayleigh number for a porous medium. It is 
worth mentioning that X > 0 corresponds to an assisting flow (heated plate), X < 0 corresponds to an 
opposing flow (cooled plate) and X - 0 corresponds to the forced convection flow. M = B 0 -J a jc p is the 
Hartmann number, Q = Q Q j '[p C p ) C ,is the heat source/sink parameter . Further, Equation 10 and 1 1 can 
be combined to give single equation: 



(1- ( p) + <p- r —f 
\.p c p) f 



f' + ff' + Qf = 0 



(13) 

Subject to the boundary condition Equation (14): 
/(0) = / D . /(«>}-> 1 
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(14) 



The physical quantity of interest is the skin friction coefficient C f , which is defined as: 

C = Zw 

where r is the skin friction or the shear stress at the surface of the plate, which is given by: 

r ~ = ^ (S) y=0 

Substituting Equation 7 into Equation 15 and 16, we obtain Equation 17: 



(15) 



(16) 



(17) 



III. RESULTS AND DISCUSION 

The nonlinear ordinary differential Equation 13 subject to the boundary conditions 14 was solved 
numerically using the shooting method. This wellknown technique is an iterative algorithm which attempts to 
identify appropriate initial conditions for a related Initial Value Problem (IVP) that provides the solution to the 
original Boundary Value Problem (BVP). The shooting method is based on MAPLE "dsolve" command and 
MAPLE implementation, "shoot" (Meade et ai, 1996). The results are given to carry out a parametric study 
showing the influences of the non-dimensional parameters, namely the mixed convection parameter X □ and the 
constant suction/injection parameter / , the Hartmann number M , the heat source/sink parameter y. 
Following Oztop and Abu-Nada (2008) we have considered the range of nanoparticles volume fraction 
(/> □ asO < <j> < 0.2 . The thermophysical properties of fluid and nanoparticles (Cu, AI2O3, TiCy used in this 
study are given in Table 1 . In order to validate the accuracy of the numerical method used, the present results for 
the skin friction coefficient /"(0) when </> = 0.1 and </> = 0.2 ,M=0.5,Q=0.5 for Cu nanoparticles and various 

values of X U are compared with those of yasin as shown in Table 2 and 3, respectively. It is clearly seen that 
the comparison shows very good agreement. Table 2 and 3 also illustrate the influence of the suction and 

injection parameter / =0.5 (suction), 0 (impermeable) and -0.5 (injection) for Cu nanoparticles and various 

values of X . The results indicate that the imposition of suction (/ > 0 ) at the surface has the tendency to 

increase the skin friction coefficient /"(0)but for the case of surface injection (/ < 0), the skin friction 

coefficient / "(0) decreases, which show a favorable agreement with the previous investigations for the case of 

impermeable surface (/ "(0)) . Figure 1 shows the variation of the skin friction coefficient (2 Pe t C . with 

X □ for different types of nanoparticles (Cu, AIO3, TiC^) when^ = 0.1. This figure shows that it is possible to 
get dual solutions of the similarity Equation 13 subjected to boundary conditions 14 for the opposing flow case 
(X < 0)with upper and lower branch solutions. Dual solutions exist fori < X > 0 , a unique solution exists 

for X = A < 0 and no solutions exist for X = X > 0 , where X is the critical value of X □ for which the 

solution exists. As in similar physical situations, we postulate that the upper branch solutions are physically 
stable and occur in practice, whilst the lower branch solutions are not physically realizable. This postulate can 
be verified by performing a stability analysis but this is beyond the scope of the present paper. On the other 

hand, it is also shown in Fig. 1 that suction (/ > 0) delays separation compared to the impermeable surface or 
injection (/ < 0) cases. This is true for all the three nanoparticles (Cu, A1203, Ti02) considered. The 
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variation of the skin friction coefficient (2 Pe ) ' C . with suction/injection parameter f B when (j> = 0.1 and 

X = -1.6 are presented in Fig. 2. This figure supports the dual nature of the solutions to the boundary value 
problem (14) and (15). For this value of A , there is a critical value / of / , at which there is a saddle-node 

bifurcation, with two solutions for (/ > / ) and no solutions for / < / < 0 . This indicates that injection 

(having (/ < 0)) limits the existence of solutions, whereas no such limit appears for suction (/ > 0), with 

both branches of solutions continuing to large values of / > 0 (suction). Based on our computations, the 

values of / 0 are f 0c = -0.515151 for A1203, f 0c = -0.49825 for Ti02 and / 0 = -0.33955 for Cu. This 

shows that the range of solutions for the injection (/ < 0) case is largest for AI2O3 nanoparticles and smallest 
for Cu nanoparticles. 




Fig. 1. Variation of the skin friction coefficient (2 p e ) C f with □ □ for different types of nanoparticles when 
□ □ = 0. 1 and various values of $ 




% 

Fig. 2. Variation of the skin friction coefficient (2 p e C f with f 0 for different types of nanoparticles when 
□ □ = 0.1 and □ D=i.6 

Figure 3-5(a-b-c-d) show the Comparison of the velocity profiles whit the effect of heat source/sink and 
magnetic field when 0<M<1, 0<Q<0.5 and <j> = 0.1 for different types of nanoparticles, namely Cu, AI2O3, 
TiC>2 respectively. The dashed line refers to the Compare article (Yasin 2013). As you see in Figs 3-(a,b,c,d) 
which refer to CU nanoparticle, the maximum point for both articles, the present work and yasin (2013) is less 
than 5 in x direction. Fig (3-c) shows that whit Q=0, the velocity profiles are in the same maximum point, it 
means that, the present results are in good agreement with those reported by yasin and we also seen that ,the 
velocity profiles decrease whit lower number of M, and the maximum point increase whit higher number of Q 
(the heat source/sink parameter). Figures 4-(a,b) which refers to AI2O3 nanoparticle have the similar behavior 
whit cu nanoparticle. 
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Fig 5(a,b,c,d) also show the results for T1O2 nanoparticle. As figures show the higher number of M (Hartman 
number) increase the velocity profiles, and the lower number of Q decrease the velocity profiles, and the suction 
parameter increases the velocity profiles and the injection parameter decreases the velocity profiles. Finally, it is 
worth mentioning that all the velocity profiles presented in Fig. 3-5 satisfy the far field boundary conditions (15) 
asymptotically. 




present work 
yasin.M H M (2013) 



Fig.(3-a) .Cu nanoparticles whit M=0.5 , Q=0.5 




Fig.(3-b). Cu nanoparticles whit M=0, Q=0. 




for M=0.5 , Q=[ 



yasin M.H.M(2013) 



Fig.(3-c) . Cu nanoparticles whit M=0.5,Q=0 



Fig. 3(a,b,c). Comparison of the Velocity profiles f\rj) for Cu nanoparticles whit yasin (2013) when 
0 = 0.1, k = -1.6 and various values of / ,M,Q 
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D 5 10 15 



Fig.(4-a) . AL 2 0 3 nanoparticles whit M=0 , Q=0.5 




Fig.(4-b) . AL 2 0 3 nanoparticles whit M=0 , Q=0.2 
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x 



Fig.(5-b) . Ti0 2 nanoparticles whit M=l , Q=0 




Fig.(5-c) . Ti0 2 nanoparticles whit M=0 , Q=0.2 




-05 1 1 1 

0 5 10 15 

Fig.(5-d) . Ti0 2 nanoparticles whit M=0 , Q=0.5 

Fig. 5(a,b,c,d). Comparison of the Velocity profiles f '(rj) for Ti02 nanoparticles whit yasin (2013) when 
0 = 0.1, A = -1.6 and various values of / ,M,Q 
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Table 1. Thermophysical properties of fluid and nanoparticles (Oztop and Abu-Nada, 2008) 



Physical properties Fluid phase (water) 


Cu 


A1203 


Ti02 






P (Kg/m3) 


999.100 


8933 


3970 


4250.000 








Table 2. Values of / '(0) for Cu nanoparticles when □ □ 


— \J. 1, ivl — V.J}^— \J.J 










Yasin M.H.M(2013) 


present 










X 


F„=0 


F 0 =-0.5 


F„=0 


F„=0.5 






-1.45 


0.39852 


1.763 


1.835 


1.796 






-1.50 


0.39263 


1.769 




1.806 






-1.516 














-1.55 


0.38391 


1.799 


1.850 


1.816 




■ 


-1.60 


0.37176 


1.789 


1.823 


1.822 






-1.65 


0.35523 










■ 


-1.70 


0.33259 


1.792 


1.843 


1.884 






-1.85 


0.30004 


1.689 










-1.90 


0.24202 


1.858 


1.856 


1.917 






-1.95 


0.18714 












-2.225 






Table 3. Values of / '(0) for Cu 


nanoparticles when □ □ 


= 0.2,M=0.5 , Q=0.5 








X 


Yasin M.H.M(2013) 
F„=0 


present 
F 0 =-0.5 


F„=0 


F„=0.5 






-1.75 


0.34746 


1.787 




1.715 






-2.00 


0.34528 


2.334 


1.738 


1.714 






-2.14460 














-2.20 


0.32409 


0.8833 


1.736 


1.711 






-2.30 


0.30311 












-2.35 


0.28845 


0.3941 










-2.40 


0.26949 














-2.50 


0.19559 












-2.60 



IV. CONCLUSION 

The aim of this present paper is to investigate numerically heat source/sink effect on the steady two- 
dimensional laminar MHD mixed convection Hiemenz flow over a vertical permeable flat plate placed in a 
porous medium and saturated by a nanofluid as considered by Ahmad and Pop (2010) and yasin (2013). We 
have extended the previous work by taking into consideration the effects of suction or injection and magnetic 
field and the effects of heat source /sink,with permeable surface. Further, the governing equations are 
transformed into ordinary differential equations and are then solved numerically using the shooting method. The 
effects of the suction or injection and the heat source/sink parameter, the mixed convection and magnetic field 
parameter and the nanoparticle volume fraction parameter on the flow and heat transfer characteristics are 
studied. In general, imposition of suction is to increase the velocity profiles and to delay the separation of 
boundary layer, while the injection parameter decreases the velocity profiles. On the other hand, the higher 
number of M (Hartman number) increase the velocity profiles, and the lower number of Q (heat source/sink 
parameter) decrease the velocity profiles and when the Q=M=0 the diagrams are the same.the range of solutions 
for the injection case is largest for A1 2 0 3 nanoparticles and smallest for Cu nanoparticles. 
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